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INTRODUCTION 


Noise production resulting from the convection of boundary-layer turbulent 
eddies over the trailing edge of an aircraft wing is regarded as an important 
source of community noise frm airframes (ref. 1 ) . An extensive experiment 
was conducted (ref. 2) to determine the physical mechanimn of this " trail ing- 
edge-noise" generation, using a two-dimensional NACA 63i-0V2 airfoil immersed 
in a uniform stream. The experimental findings indicate that the flow struc- 
ture of the boundary layer and near wake in the vicinity of the airfoil trail- 
ing edge play a central role in noise generation. The peak of the broadband 
trailing-edge noise that is generated scales with convection speed and boundary- 
layer thickness. Evidence shows that the intensity of the noise radiated 
depends on the changes in turbulence and mean-flcw properties which occur at 
the trailing edge. Therefore, there is a need to understand the structure of 
mean and turbulent flow near the trailing edge in order to shed more light on 
the trailing-edge-noise production mechanism. 

Prom the aerodynamic point of view, flow near the trailing edge of an air- 
foil is an important aspect in airfoil design. The trailing-edge flow controls 
the circulation around the airfoil and, hence, influences the entire flow field. 
Although detailed flow measurements at low and moderate free-stream Mach numbers 
are available for the near wake of flat plates (ref. 3) and airfoil cascades 
(ref. 4) , the same is not true for an isolated airfoil. Detailed turbulent- 
flow measurements in the boundary layer and near wake of an isolated airfoil 
at low and moderate free-stream Mach numbers have received proper attention only 
recently (ref. 5) . 

In support of the present study of airfoil trailing-edge noise, a computa- 
tional analysis was initiated (ref. 6) to predict the flow-field details near 
the airfoil trailing edge. The analysis is based on a finite- element formulism. 
A parabolized form of the Navier-Stokes equations is solved in conjunction with 
analysis of a two-dimensional, inviscid potential flow. The computation code 
is initialized with distributions of mean flow and Reynolds stress measured 
upstream (90-percent chord) of the trailing edge. Detailed flow-field predic- 
tions were generated for the NACA 63i-012 airfoil at zero angle of attack. 

The purpose of this paper is to present the flow-field investigation 
obtained with the two-dimensional NACA 63^ -01 2 airfoil. Measurements made 
Include profiles of mean streamwise velocity and distribution of Reynolds 
stress in a region of 20-percent chord centered at the airfoil trailing edge. 
Both naturally developed and tripped boundary- layer flows were investigated. 

The Reynolds number based on airfoil chord was 1.25 x 10®, and the angle of 
attack was set at zero. These measurements provide the required flow-field 
details for understanding the trailing-edge-noise production process and also 
serve As the data base for validation of the computation code. Comparisons 
between measured and ccanputed results for the naturally transitioned flow were 
previously reported in reference 6. The computed results agreed with the mea- 
surements reasonably well. 



Identification of commercial products in this report is used to adequately 
describe the model. The identification of these commercial products does not 
constitute official endorsement, expressed or implied, of such products or manu- 
facturers by the National Aeronautics and Space Administration. 


SyMBOLS AND ABBREVIATIONS 

airfoil chord, m 

local skin-friction coefficient 

diameter of hot-wire sensor, m 

shape factor, 6*/6 

length of hot-wire sensor, m 

tur bulent- kinetic-energy flux, (u')^ + (v')^» mVs^ 
total turbulent-kinetic-energy flux, (u')^ + (v')^ + (w')^f 
Reynolds number based on boundary-layer thickness 
measurement station 

streamwise (tangential) mean velocity, m/s 
free-stream velocity, fly's 

streamwise, turbulent, root-mean- square velocity, m/s 

Reynolds shear stress. Pa 

frictional velocity, m/s 

normalized wall velocity, U/u-^, m/s 

transverse, turbulent, root-mean-square velocity, m/s 


wake function, 2 sin^ 



w' spanwise, turbulent, root-mean-square velocity, m/s 

X streamwise distance along chordline, m 

y transverse distance normal to airfoil surface (or normal to wake 

centerline) , m 
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UxY 

y"'' normalized transverse distance from wall in wall coordinate, — 

V 

y* normalized transverse distance from wall, y/6 

a wake coefficient 

6 boundary-layer thickness, m 

6* boundary-layer displacement thickness, m 

6 angle between mean flow and normal vector of wire, deg; and momentum 

thickness, m 

V kinematic viscosity, m^/s 

1 I absolute value 

Subscript: 

max maximum 

min minimum 


EXPERIMENTAL APPARATUS, INSTRUMENTATION, AND PROCEDURES 

The experiment was performed in the Langley Aircraft Noise Reduction Labor- 
atory. The facility is essentially a continuous-flow, open-section, open-circuit 
wind tunnel in an anechoic chamber. The rectangular nozzle used to generate the 
flow had exit dimensions of 0.30 m x 0.46 m and an overall, area contraction 
ratio of 34. Calibration measurements showed that the exit mean flow was uni- 
form within 1.0 percent of the nozzle exit. The streamwise turbulent intensity 
was less than 0.5 percent at the center of the nozzle exit. Details of the air 
supply system and the anechoic chamber are available in reference 7. 

Airfoil and Test Section 

An NACA 63i -01 2 symmetrical airfoil with 0.61 -m chord and 0.30-m span were 
used in the experiment. The airfoil was machined from stock aluminum blocks 
and the cusped trailing edge was made knife sharp. This airfoil is one of two 
identical airfoils used in a previous experiment (ref. 2) and was instrumented 
with a number of flush-mounted transducers to measure fluctuating surface 
pressure. 

The test section, which provided the two-dimensional flow for the airfoil, 
consisted of two parallel sideplates mounted along the short edges of the nozzle. 
The airfoil was installed between the sideplates with the leading edge at 
25-percent chord downstream of the nozzle exit. A schematic of the test setup 
is given in figure 1 (a) . The corresponding photograph is shown in figure 1 (b) . 
Other details of the test section may be found in reference 2. 
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Artificial boundary tripping was achieved by gluing a 0. 076-mm-diaineter 
Steel wire across the entire span of the airfoil at 6-percent chord f ran the 
leading edge. To maintain the symmetry of the flow, identical tripping arrange- 
ments were used on both surfaces of the airfoil. The purpose of the tripping 
was to artificially induce boundary-layer transition near the leading edge so 
that a thicker boundary layer resulted at the trailing edge. The thickened 
boundary layer was expected to simulate the higher Reynolds flow conditions 
prevailing over a full-scale aircraft wing. 


Hot-Wire Probe and Instrumentation 

At the test Reynolds number of the airfoil, the boundary layer was esti- 
mated to be of the order of 10 mm. In order to obtain flow details in the wall 
region, it was necessary to employ a very small x-wire probe so that adequate 
spatial resolution and close approach to the wall were possible. A special, 
subminiature x-wire probe was developed for this study. The tungsten wire used 
had a 2.5-um diameter, and its length to diameter ratio il/d was about 200. 

The lateral spacing between the two conponent wires was 0.50 mm. A schematic 
drawing and a photograph of the x-wire probe are shown in figure 2. The small 
metallic pin mounted on the back of the probe body was used as a contact indi- 
cator to provide an accurate initial position of the probe in the boundary- layer 
measurement. The offset distance between the contact point and probe prong, 
typically about 0.02 mm, was measured with a high-power microscope. Most mea- 
surements were made with the special x-wire probe. Some measurements, however, 
were made with standard DISA^ 55P63 x-wire and 55P1 4 single-wire probes. Both 
probes use 5-ym platinum-plated tungsten wire with Jl/d = 250. Measurements 
made by the three probes agreed within 5 percent for mean velocity and Reynolds 
stress. Two channels of hot-wire instrumentation were used for the measurements. 
Each channel consists of a constant-temperature hot-wire anemoneter with a 5:1 
bridge, a polynomial linearizer, and a filter. 

The hot-wire probe was traversed in a direction normal to the airfoil sur- 
face at different chordwise stations. Most of the traverses were made in the 
midspan plane of the airfoil, although some traverses were also made at differ^ 
ent spanwise stations to check the two-dimensionality of the flow. The major 
component of the traverse rig was a remotely controlled traversing mechanism 
driven by a stepper motor and a sweep drive unit. The resolution of the tra- 
verse was 0.02 mm. During the measurement, the traverse was moved in discrete 
steps. The size of the steps was varied depending on the probe location, fine 
step size being used in the flow region near the wall and coarse steps in the 
outer region of the boundary layer. Typically, 1 00 data points were obtained 
in a flow region of one boundary-layer thickness. 

Hot-wire probes were calibrated with a free-jet calibration apparatus. The 
overheat ratio used was about 0.8. For x-wire probes, each component wire was 
calibrated with the wire normal to the flow. Typically, 12 calibration points 
were sufficient to accurately define the wire response within the measurement 
range (0 to 35 m/s). The calibration data were then curve fitted to a fourth 


^ Disa Electronics, division of Disamatic, Inc. 
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degree polynomial by the method of least squares. The coefficients of the 
polynomial were used to set the linear izer. The inclination angle of the wire 
with respect to the true probe axis was also calibrated and was found to be 
within 2° of the desired 45° angle for all x-wire probes. A typical angular 
response of a component wire is shown in figure 3. The measured response is 
compared with the law of cosines which is normally assumed for the angular 
response of an inclined wire. It is seen that the measured response closely 
follows the law of cosines. The variation from the law of cosines is within 
2 percent for 6 between 0° and 70°. This means that, if the x-wire is aligned 
symmetrically to the mean flow, a maximum measurement error of 4 percent could 
result if the instantaneous- flow incidence angle is within ±25° of the true 
probe axis. This maximum error would pccur if the root-mean-square flow fluc- 
tuation is about 45 percent of the mean velocity. Based on Klebanoff's measure- 
ment (ref. 8) in the wall region of the turbulent boundary layer over a smooth 
flat plate, flow fluctuations of this magnitude occur only when y'*' = u^y/v is 
less than 30. For the present study, the closest approach to the wall for 
x-probes was about y"^ = 40, Therefore, the departure of the wire response 
from the law of cosines was not considered significant. By use of the x-wire 
probe calibration procedure just discussed, the velocity measured by the x-wire 
and single-wire probes is canpared with that obtained with a pitot tube on the 
centerline of the calibrator nozzle. (See fig. 4.)- Over most of the measure- 
ment range, x-wire results are within 5 percent of the pitot-tube measurements. 
For the single wire, the agreement was within 1 percent. 

One of the major problems in the interpretation of hot-wire signals is 
that they are sensitive to variations of mean-flow temperature. In the present 
experiment, the actual drift of mean-flow temperature during the measurement 
was small (within 1°C) . The problem encountered was that the temperature of 
the flow in which the probe was calibrated differed from that of the flow being 
measured. Since probes were always calibrated just prior to the measurement, 
this temperature difference never exceeds 5°C. Therefore, a correction method 
developed by Bearman (ref. 9) , suitable for correcting bridge voltage with 
small, mean- flow temperature changes, was used. This, correction method was 
checked with the calibration data obtained at several different mean-flow tem- 
peratures. The accuracy of the method was found to be within 1 percent, a 
conclusion that was also reached by Bearman. If mean-flow temperature differs 
during the measurement from that used at calibration, Bearman' s method is used 
to compute the full-scale flow voltage at 9 - 0°. The bridge voltage at maxi- 
mum mean velocity for 6 = 0° is required in order to span the linearizer. 


Data Reduction 

Mean velocities and Reynolds stresses of the flow were deduced frcxn suit- 
ably amplified and linearized signals obtained frcm two hot-wire channels via 
a digital computer. Instantaneous signals were first digitized at a rate of 
TOO kHz, and then the axial and transverse components were obtained by adding 
and subtracting the two signals, respectively, at each probe location by assum- 
ing a cosine-law angular response. The mean and mean-square flow quantities 
were computed by coatinuous ensemble averaging over the data blocks. The data 
acquisition at each probe position terminates when both the mean and mean-scjaare 
values converged within specified tolerances or when the maximum allowed number 
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of averages was exceeded. In the latter situation, the probe position was 
flagged and the deviation from the mean recorded. The measurement was repeated 
when more than 1 0 percent of the data points were flagged. The tolerance for 
convergence for mean velocity was 1 percent and for root-mean-square turbulence 
quantities was 2 percent. For mean velocity measurements, ^11 data satisfies 
the convergence criteria. For fluctuating quantities in a typical probe tra- 
verse, convergence was usually met for 95 percent of the data points. 


Experimental Procedures 

For each test run, the tunnel flow was first stabilized at 30.5 m/s. The 
hot-wire probe was then moved to the uniform-flow region away fran the airfoil. 
With the probe properly aligned to the flow, the cold resistance was measured. 
The flow temperature deduced frcan measured cold resistance was used to ccmpute 
bridge voltage corresponding to maximum flow velocity. For a single-wire probe, 
the same voltage at 0 = 0° could be measured directly. A ccmparison between 
the computed voltage and the directly measured value thus served as an indepen- 
dent check for the probe alignment. With an x-wire probe, it was not convenient 
to align each wire normal (0 = 0°) to the test flow. For the procedure used, 
the bridge output voltages were measured with the probe aligned to the flow 
(0 = 45°). The computed voltages at 0 = 45° were compared to the measured 
values as a check of probe alignment. The required, full-scale, flow-bridge 
output voltage at 0 = 0° can be obtained either from Bearman's method or from 
the directional response of the wire. The two usually agreed within 2 percent. 

For each boundary-layer traverse, the probe axis was aligned tangential 
to the airfoil surface at the initial position of the traverse by use of a 
high-power telescope with a resolving power of 0.01 ram. The precise location 
of the probe, however, was provided by the contact pin arrangement discussed 
earlier. The probe was traversed in a direction normal to the airfoil surface 
in the raidspan plane of the airfoil. For wake measurements, the probe was 
aligned tangential to the chordline of the airfoil. 


RESULTS AND DISCUSSION 
Verification of Measurement Technique 

Since x-wire measurements are subject to more sources of experimental error 
(such as probe orientation, probe geometry, wall proximity interference, and 
mismatching of two component wires) than the single-wire measurements, a com- 
parison between corresponding measurements obtained by the two different types 
of probes is indicative of the reliability of the present measurement. Shown 
in figure 5 are comparisons of normalized, streamwise, mean velocity U/U^^^ for 
the naturally developed boundary layer and for the artifically tripped boundary 
layer. It is seen that the agreement is very good, within 5 percent. The com- 
parison shown in figure 6 is of normalized, streamwise, root-mean-square (rms) 
turbulent velocity u’/Uo measured fran the x-probe and the rms fluctuating 
velocity measured by the single-wire probe. It should be noted that, for a 
single wire which responds to velocity normal to the wire, the measured results 
are expected to be higher than the true u' because of strong lateral and 
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transverse velocity fluctuations. Bradshaw (ref. 10) discussed in detail the 
single-wire response in a turbulent boundary layer. For a turbulent boundary 
layer over a smooth flat plater he estimated the maximum error for single-wire 
measurement of u' to be 2 percent. The actual error depends on the probe 
location in the boundary layer and the turbulence structure of the flow. Frem 
figure 6, it is seen that the single-wire result is about 7 percent higher than 
the x-wire result, suggesting a 5-percent measurement error with the x-wire 
probe. Hence, the maximum measurement error on Reynolds shear stress -u'v* 
could approach 10 percent. 

The symmetry of the test flow was deduced from flow measur emen ts made in 
the near wake of the airfoil. Shown in figure 7 are U and -u'v' distribu- 
tions across the airfoil near the wake at x/c = 1.0029. The symmetry is within 
2 percent for both U and -u'v* . Selected traverses made across the entire 
wake flow at downstream stations indicated similar results. Thus, most of the 
wake traverses were made on the same side where the boundary- layer measurements 
were performed. 


Determination of Mean-Flow-Integral Properties 

To determine the boundary-layer mean-flow-integral properties Cf and 6, 
the measured velocity profiles were curve fitted (least-square method) to the 
wall-wake law as proposed by Coles (ref. 11): 


u+ = 


a 


In y+ + 5 + — — w(y*) 

0.41 0.41 


(1) 


where the wake function (law of the wake) w(y*) is given by 


w(y*) = 2 sin^ f- y*j 


( 2 ) 


The three parameters ux, a, and 6 are determined frcxn the fitting procedure. 
Shown in figure 8 are typical comparisons between the measurement and Coles' 
law of the wall. 


u+ = In y+ + 5 

0.41 


(3) 


It is seen that for both tripped and untripped boundary layers, the agreement 
is very good beyond transition in the near-wall region (30 < y"*" < 300). In the 
outer flow region,, the measurement deviates frcan the wall law,, more so for the 
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tripped than the untripped boundary layer. This indicates that artificial 
tripping increases the wake component in the boundary layer. The corresponding 
comparison with the wake law is given in figure 9. Agreement between the data 
and Coles' wake function is fair for the tripped boundary layer but rather poor 
for the untripped boundary layer, presumably because of the lack of wake con- 
ponents in the untripped boundary layer. 

The local skin-friction coefficient Cf was computed fron u^ and by 


Cf 



(4) 


The momentum thickness 6 and displacement thickness 6* for the wake 
flow were obtained by integrating the measured velocity profiles in y using 
a Runge-Kutta third-order numerical scheme. Since measurements cannot be 
obtained for y"*" < 40, direct integration over only the measured velocity can 
introduce large errors in 0 and 6* for boundary- layer flow. This is because 
the velocity distribution near the wall (y"*" < 40) has the greatest eontrlbution 
to the integrals of 9 and 6*. in view of the good fit for the present mea- 
surements obtained with Coles' law of the wall, the mean velocity was extrap^ 
dated to the wall through Coles' law (eq. (3)) and the linear law given by 

u"*" = y"*" (5) 


The two flow regions intersect at yj - VO. 8. The velocity extrapolation was 
done with equation (5) from the wall to yj and with equation (3) from yj 
to y+ - 40. Although the transitional behavior fron linear region to log 
region is not taken into account in the extrapolation, it is felt that the 9 
and 6* obtained with this extrapolation should be much closer to the true val- 
ues than otherwise. The improvements obtained with this extrapolation for the 
untripped flow were about 1 5 to 20 percent in 9 and 25 to 30 percent in 5 . 
For the tripped flow, the improvements were about half of those given above. 

Measured flow properties for both naturally developed flow and arti- 
ficially tripped flow are presented in the following sections. Mean veloc- 
ity and rms turbulent velocities u^ and v' are normalized to free-stream 

velocity. Reynolds shear stress -u'v' and turbulent kinetic energy flux 
q2 are normalized to the square of free-strea m ve locity. Also presented are 
the turbulent-shear correlation coefficient -u'v'/u'v' and the ratio of 
Reynolds shear stress to turbulent kinetic-energy flux lu'v'|/q^. 


Naturally Developed Flow 

The development of streamwise mean velocity for the untripped flow is given 
in figure 10, over the flow region 0.9 5 x/c = 1.1. The chordwise locations 
of the measurement stations are given in table 1 . For boundary- layer flow 
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(x/c < 1.0), the mean velocity profiles are similar. By proper integration 
and/or fitting with Coles' wall-wake law, integral properties 6*, 0, H, 

and Cf were obtained from these mean velocity profiles. Results are given 
in figure 11 . Both 6* and 0 are essentially constant at 1 mm and 
1.4 mm, respectively, in the boundary layer, with a slight initial increase at 
x/c ® 0.94. In the near wake, both quantities exhibit maxima of about 1.1 mm 
and 1.5 mm, respectively, at x/c ^ 1.005. The shape factor H is nearly con- 
stant (®1.39) in the boundary layer. Downstream of the trailing edge, H 
is seen to exhibit an initial sudden reduction (to about 1.3), followed by a 
gradual decrease towards unity as expected for the far wake. The local skin- 
friction coefficient Cf exhibits some scatter. The general trend is a slight 
reduction of skin friction as the airfoil trailing edge is approached. The 
average value of Cf is about 0.0039 over the last 10 percent of the airfoil 
chord. 


Distributions of u', v', -u'v', and q- are presented in figures 12 

to 15. The peak, streamwise, turbulent intensity u'/U^^ is about 7.5 percent 
in the boundary layer and increases to about 8.0 percent just downstream of the 
trailing edge. This is followed by a gradual decay farther downstream. Also 
noted is a spike in u'/U^^ in the slipstream where the two boundary layers from 
opposite sides of the airfoil meet. This spike migrates outward and merges with 
the broader peak observed in the boundary- layer region. The position of the 
broad peak in u' is found to progressively move away fron the chordline 
indicative of the spreading of the mixing region in the wake, as shown in fig- 
ures 12(a) and 12(b). The transverse turbulent intensity v'/U^ has a peak 
value of 6 percent at x/c = 0.9. The peak value occurs very close to the sur- 
face and decays continuously with downstream location in the wake. Of particu- 
lar interest is the reduction in v'/U^g across the trailing edge, (Ccxnpare 
sta 11 with sta 12 in fig. 13.) At x/c - 1.1, the peak value is reduced to 

about 3.5 percent. The distribution for v' exhibits a spike in the slip 

stream similar to the u' distribution. This spike, however, remains on the 
chordline of the airfoil and decays gradually with downstream distance. Hence, 
the peak value of near wake remains on the chordline. 

Distribution of the Reynolds shear stress is given in figure 14. At 
x/c = 0.9, the normalized, peak shear stress is about 0.0012. The peak value 
then increases rapidly to 0.0017 at x/c - 0.9583, followed by a decay to its 
initial value at x/c = 0.9896. Prom x/c = 0.9896, the peak value of the 

normalized shear stress decreases. In the slipstream, a spike is also observed. 

The spike diffuses outward and eventually merges with the broad peak in the 
outer region into a smooth distribution. The diffusion of the slipstream spike 
for -u'v' is seen to be similar to that for u'. Fran boundary layer to the 
near wake, the Reynolds shear stress shows a reduction which occurs immediately 
downstream of the trailing edge. 


Distributions of turbulent-kinetic-energy flux q^ are shown in figure 15. 
It should be noted that the q^ defined here is not the total turbulent- kinetic- 
energy flux, because the spanwise velocity fluctuation w' was not measured. 

An estimate of the total turbulent-kinetic-energy flux (qt^ - u'^ + v'^ + w'^) 


may be obtained if it is assumed that w'^ - 0.5(u'2 + v'^), as was done in 
reference 5. In that case, q^^ = 1.5q2 and the total turbulent-kinetic- 


9 



energy flux would be 50 percent higher than the value of plotted. 

As shown in figure 1 5f remains nearly constant for the four initial mea- 
surement stations (0.9 ^ x/c ^ 0.9688) with ® 0.09. From x/e = 0.9688 


to x/c as 1 , a gradual decrease in q^ is noted with qVu^^ ® 0.07 at the 
trailing edge. For all boundary-layer stations, q2 peaks very close to the 
airfoil surface. In the near-wake region, note that a high q2 concentration 
occurs in the slipstream. The effects of the lateral turbulence mixing in 
the wake region, however, diffuses this high concentration of turbulent kinetic 
energy over a distance of about 6-percent chord. At x/c = 1.1, the location 
of maximum q^ is seen to be displaced to the middle of the wake mixing layer. 


The distribution of turbulent- shear correlation coefficient is shown in 
f igur e 16. This is obtained by dividing the measured Reynolds shear stress 
-u'v' by the product of the root-mean-square values of u' and v' . The cor- 
relation coefficient distribution measured in the boundary-layer flow typically 
has a positive maximum near the airfoil surface and remains nearly constant at 
the maximum value for the inner part of the boundary layer. As the outer edge 
of the boundary layer is approached, the correlation coefficient drops rapidly 
to zero. If the boundary-layer flow is generated by a very low, turbulent, uni- 
form free-stream flow, the correlation coefficient should remain nearly zero 
outside the boundary layer. In the present experiment, the free-stream flow 
was provided by the potential core of a free jet. At the airfoil trailing edge, 
the thickness of the free-jet shear layer was about 7.6 cm and the distance from 
the airfoil trailing edge to the inner edge of the free-jet shear layer is about 
11.4 cm. Even thou gh the turbulence intensities u' and v' and Reynolds 
shear stress -u'v' measured in the free-stream of the airfoil are extremely 
low (see figs. 12 to 14), the distribution of correlation coefficient in fig- 
ure 16 reveals the presence of the free-jet shear layer. This is indicated by 
the reversal of the sign of the correlation coefficient as the measurement is 
made farther away fron the outer edge of the boundary layer in the free-stream 
flow. 


Similar observations can also be made in the near-wake region of the air- 
foil. It is pertinent to note at this point that, although the distribution 
of the correlation coefficient rev eals the presence of the free-jet shear layer, 
the magnitudes of u', v', and -u'v' are extremely low and would not affect 

the validity of the same measurements made inside the boundary layer, fhe 
maximum value of the correlation coefficient varies from 0.25 at x/c = 0.9 
to 0.4 at the trailing edge. In the near-wake region, a further increase is 
noted and a value of about 0.5 for the correlation coefficient is reached at 
x/c — 1.1. Comparing the shear correlation across the trailing edge (sta 11 
and sta 12 in fig. 16), it is seen that the shear-correlation profile displaces 
away from the chordline in the near wake. Depending on the Reynolds number 
based on 6, correlation coefficients of about 0.5 have been observed in a fully 
developed, turbulent boundary layer over a flat plate under zero pressure gradi- 
ent (ref. 8) and in a fully developed, turbulent, two-dimensional channel flow 
(ref. 12). The fact that there exists a region near the wall over which the 
correlation coefficient remains nearly constant is indicative of the presence 
of a constant stress layer as expected frcsn theoretical arguments (ref. 13). 

In the near-wake region, the extent of the constant-stress region is seen to 
broaden with downstream distance and is further displaced frcxti the chordline 
of the airfoil. 
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Figure 17 is an illustration of the distribution of |u'v'|/q2. This ratio 
is a measure of the relative importance of turbulent-kinetic-energy production 
in the flow and the turbulent-kinetic-energy content in the flow. For most 
turbulent boundary layers, the ratio is nearly a constant in the outer region 
of the boundary layer. For example, the ratio is about 0.167 for a turbulent 
boundary layer over a flat plate under zero pressure gradient ( ref. 1 4) . For 
the present airfoil, there are indications in figure 17 that |u'v'|/q2 in the 
boundary layer takes on a value of 0.13 at x/c = 0.9, increases to a maximum 
of 0.20 at x/c = 0.9844, and then decreases to 0.18 at the trailing edge. In 

the ne ar wa ke, an initial decrease is noted just downstream of the trai ling edge 

with [u'v'l/q^ = 0.18, followed by a rather gradual increase to lu'v'|/q^ = 0.2 
at x/c = 1 .1 . Since the value of q^ used is just the sum of mean squares 
of the streamwise and transverse velocity fluctuatio ns, i t represents only 2/3 
of the total turbulent kinetic energy qt^. Thus, lu’v'l/q^^ would be about 
0.13, a value somewhat lower than that established for a flat-plate turbulent 
boundary layer. 

It is also interesting to note that the region over which |u’vVl/q2 is 
constant exists only over a limited extent in the boundary layer. The sharp 
drop-off of this ratio in the outer region of the boundary layer is also not 
o bser ved for a flat-plate turbulent boundary layer. In the near wake, however, 
|u'v’|/q2 is nearly constant over most of the outer region of the wake and 

decre ases sharply to zero near the centerline of the wake. The nonzero value 

of lu'v'l/q^ observed in the free-stream flow is due to the influence of the 
f ree- jet shear layer as discussed earlier. It should be noted, however, that 
|u'v'| and q^ are extremely low in the free stream. 

In regard to the dependence of tr ailing-edge-noise production on the 
boundary-layer/near-wake turbulence properties near the edge, it has been 
observed for the untripped flow tha t u' increases just downstream of the 
trailing edge while v' and -u'v' decrease and the shear correlation profile 
displaces away from the chordline. These observations support the contention 
made in reference 2 that the coherent eddies undergo an adjustment across the 
trailing edge in the process of producing trailing-edge noise. 


Artificially Tripped Flow 

In the second part of the investigation, the boundary-layer growth along 
the airfoil was enhanced by artificial tripping. Distributions of streamwise 
mean velocity are given in figure 18. It is seen that the boundary layer is 
thickened by the leading-edge tripping (compare fig. 18 with fig. 10). The 

profiles are nearly identical in the boundary layer and also in the outer 
part of the wake. This is indicative of the establishment of a fully developed, 
equilibrium, boundary-layer flow near the airfoil trailing edge. The integral 
mean-flow properties deduced fran the U distributions are presented in fig- 
ure 19. At the first measurement station (x/c = 0.9) in the boundary layer, 
the values of 6 and 6* are 1 .87 mm and 2.95 mm, respectively. The boundary- 
layer displacement thickness 6 exhibits an initial increase, reaching a maxi- 
mum value of 3.05 mm at x/c - 0.98, followed by a reduction towards the trail- 
ing edge. Similar behavior is noted for 8. Both 6 and 6* are seen to 
increase to a maximum value just downstream of the trailing edge at x/c = 1 .005, 
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similar to the untripped flow. This behavior is believed to be typical of flow 
transition from bounded, boundary-layer shear flow to free-shear, near-wake flow. 
The shape factor was reasonably constant at 1 .55 at initial measurement locations 
fr<xn x/c =0.9 to x/c = 0,98. Close to the trailing edge, H starts to 
decrease, reaching a value of 1 .51 at x/c ^ 0.995. This decrease continues 
into the near wake and tends to approach unity asymptotically, as expected, for 
the far wake. The major difference noted in the mean-flow integral properties 
between the tripped and untripped boundary layers lies in the local skin- 
friction coefficient. While a slight decrease of Cf is noted towards the 
trailing edge of the untripped flow, the tripped boundary layer Indicates a 
definitive opposite trend. The averaged values of Cf over the last 10-percent 
of airfoil chord are 0.0025 and 0,0039 for the tripped and untripped boundary 
layer, respectively. 

The distributions of and v'/U^ are given in figures 20 and 21 , 

respectively. The u'/U,,;, distributions measured in the boundary layer are 
similar. The common peak value is about 7 percent, which occurs at y/6* = 1,3, 
In the near-wake region, the outer part of the u’/Ugg distributions are also 
nearly identical to that measured in the boundary layer. In the slipstream 
downstream of the trailing edge, a hump appears, as previously occurred for the 
untripped flow. The lateral diffusion of the hump, however, is more rapid than 
for the untripped case. Measured v'/U^ distributions show trends similar 
to those for the u'/U„ variation. The peak value in the boundary layer and 
the outer part of the wake is about 5.5 percent. The general behavior of the 
slipstream spike in the wake region is similar to that noted for the untripped 
flow. Contrary to the untripped flow, there is no change observed in u' or 
v' across the trailing edge other than the spikes in the slipstream. 

Distributions of Reynolds shear stress and turbulent kinetic energy for the 
tripped flow are presented in figures 22 and 23, respectively. Similar distri- 
butions of -u'v'/Uoo^ are noted in the boundary-layer region with a peak value 
of 0.0014. A slight increase in the peak value is found in the initial part 
of the measurement r egio n (x/c ® 0.9) . Downstream of the airfoil trailing edge, 
the peak value of -u’ v'/Ooo^ increases to 0.0017 and remains at this value to 
x/c = 1.02 followed by a decay to 0.001 2 at x/c = 1 .1 . The variation of 

is similar to that of the Reynolds shear stress. The peak value is about 
0.008 which occurs at 0,36 from the airfoil surface. In comparison with cor- 
responding measurements made for the untripped flow, it is seen that the peak 
values are comparable for the two flows, but the tripping displaces the region 
of maximum, turbulent kinetic energy farther away from the surface. In the 
near -wake region, the distribution in the outer part of the wake is 

again similar to those measured in the boundary layer. A decrease of turbulent 
kinetic energy is seen to occur for x/c > 1 .06. The high level of q^ in the 
slipstream is seen to diffuse more rapidly for the tripped flow than that for 
the untripped flow. 

Shown in figure 24 is the distribution of turbulent-shear correlation 
coefficient and in figure 25 is the ratio of turbulent energy production and 
turbulent kinetic energy. The turbulent shear correlation remains nearly con- 
stant in the boundary-layer region. The maximum value of -u'v'/u'V is about 
0.32. A slight increase in turbulent shear correlation may be noted from 
x/c =0.9 to x/c = 0.9417. In the wake region, the maximum value of the 
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correlation coefficient is about 0.38, and a slight decrease is noted at 
x/c = 1 .1 . In comparison with the corresponding measurements made in the 
untripped flow, it is seen that the effect of tripping is the elimination of 
most of the streamwise variations of the shear correlation coefficient. More- 
over, -u'vVu'v' is seen to maintain a constant value over most parts of the 
boundary layer and/or wake for the tripped flow. On the other h^nd, tripping 
reduces the mag nitu de of the correlation coefficient. The ratio of Reynolds 
shear stress [u'v'l and the turbulent kinetic energy (see fig. 24) is 

seen to follow a behav ior similar to that noted for the shear correlation coef- 
ficient. The ratio | u ' v ' | /q^ remains nearly constant over most parts of the 
boundary layer and/or wake with very little streamwise variation. In the bound- 
ary layer, the ratio is about 0.16, and in the wake, it is about 0.18. Both 
values are lower than those found for the naturally develope d flo w. The influ- 
ence of the free- jet shear layer can be inferred from the -u'v'/u'v' and 
lu'v'|/q^ distributions outside th e bo undary layer and wake flow. As pointed 
out previously, the magnitudes of u'v', u', v', and q^ are extremely small 

outside the boundary layer and wake; hence, the existence of the free-jet shear 
layer does not influence the accuracy of the measurements made in the boundary- 
layer and wake flows. 

Of particular significance, the changes in turbulence properties across 
the trailing edge for the tripped flow exhibit a different behavior from that 
observed for the untripped flow. The only change noted for the tripped flow is 
in the Reynolds shear stress which showed an increase downstream of the trailing 
edge. As discussed in the previous section, Reynolds shear stress and v', for 
the untripped flow, were found to decrease in the immediate downstream region 
of the trailing edge, while u' showed an increase there. These contrasting 
behaviors are believed to be a result of increased fine-scale turbulence due 
to tripping as well as a result of loss of correlation of the coherent eddies 
in the boundary layer. As far as the trailing-edge-noise mechanism is con- 
cerned, more refined measurement such as a conditioning sampling technique may 
be necessary in order to resolve the difference observed between the two types 
of flow. 


Comparison of Present Untripped Results With Measurements of 
Turbulent Boundary Layer Over a Flat Plate 

The flow field around the present airfoil is affected by the streamline 
curvature and variable pressure gradient. The airfoil surface near the trailing 
edge (x/c > 0.9) is concave, and a mild positive pressure gradient prevails near 
the trailing-edge region. It is of interest to qualitatively compare the present 
results for the naturally developed flow with those obtained for turbulent flow 
over a flat plate with zero pressure gradient. In this way, the influences of 
mean-flow pressure gradient and streamline curvature on the development of a 
boundary layer may be delineated. 

The mean, streamwise velocity distribution for the present airfoil is seen 
to obey the law of the wall for y"*" ^ 300. (See fig. 8.) In the outer region, 
a weak wake component is also noted. For a flat-plate boundary layer, the law 
of the wall generally holds (ref. 13) for y+ up to 1000. The measured, local 
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skin-friction coefficient is about 0.0039 for the untripped airfoil as compared 
to a value of 0.0036 predicted by the empirical formula of Ludweig and Tillmann 
(see Hinze, ref. 1 3) for the flat-plate boundary layer 


Cf 


0.246 • 10-0 • 6 . 



(6) 


A 

valid for — from 1 O'* to 10^. The shape factor was deduced from the data to 
V 

be 1.39 for the untripped-airfoil boundary layer/ compared to a value of 1.39 
predicted by Hamma's empirical formula (see Hinze/ ref. 13) for the flat-plate 
boundary layer 


H = 


- 6.3 



(7) 


For the turbulent quantities/ the present measurements can be compared with 
those reported by Klebanoff (ref. 8) for a fully developed, turbulent boundary 
layer over a flat plate. The Reynolds number based on boundary-layer thickness 
6 was 7.5 X 10^ for Klebanoff 's data and is 2 x l o'* at the trailing edge of 
the present, naturally developed boundary layer. The value of u'/U^q at 0.1 6 
was about 7.3 percent as determined by Klebanoff s' measurement, compared to a 
value of 7 percent obtained at the airfoil trailing edge. The ratio v'/U 
reached a peak of about 4 percent for the flat-plate boundary layer. The air- 
foil measurements indicate a 5-p erce nt ^ak near the trailing edge. The nor- 
malized Reynol ds s hear stress -u'v'/U„ and turbulent-shear correlation 
coefficient -u'v'/u'v' were 0.014 and 0.5, respectively, from Klebanoff 's mea- 
surements^ For the airfoil, the corresponding values are 0.017 and 0.4. The 

ratio lu'v'|/qt^ exhibited a flat peak of 0.166 for the flat plate and 0.13 
for the airfoil. 

In general, the boundary-layer measurements made near the trailing edge 
of the untripped airfoil follow closely the trends observed in a flat-plate 
boundary layer under zero pressure gradient. This observation is indicative 
of the weak influences of the streamline curvature and the adverse pressure 
gradient on the airfoil boundary layer. It should be pointed out, however, 
that the observed variation of flow properties in the flow direction as dis- 
cussed in Naturally Developed Flow would be absent for a fully developed, tur- 
bulent boundary layer over a flat plate. 


Effect of Tripping on Airfoil Boundary Layer 

As was stated previously, the purpose of boundary-layer tripping is to 
enhance the flow development and to promote the mixing so that a thicker bound- 
ary layer is obtained near the trailing edge. It should be pointed out that 
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a tripping device such as the one used in the present study, which is common 
to this type of study, may alter the turbulence structure downstream, as cau- 
tioned by Willmarth (ref. 1 4) . This is especially true for large-scale distur- 
bances in the flow. Nevertheless, the thicker boundary layer resulting from 
the tripping does facilitate measurement and does increase experimental accu- 
racy. In addition, the transition is artificially fixed by the tripping wire, 
which causes minimization of the nonuniformities in the spanwise direction of 
the flow. The measurements obtained for the tripped flow are considered repre- 
sentative of those which would exist if an airfoil of longer chord were used. 

The most fundamental difference found between the naturally developed flow 
and the tripped flw is the lack of streamwise variation of the latter. The 
artificially tripped flow exhibits a high degree of uniformity in all normalized 
flow properties, both in mean flow and in turbulence quantities. For the mean- 
flow field, tripping increases the wake canponent in the outer region, reduces 
the local skin-friction coefficient, and slightly increases the shape factor. 

The value of R§ for the tripped flow is about 5.1 x 10^, a factor of 2.5 
increase over the untripped flow. Laufer (ref. 1 2) reported turbulence measure- 
ments in a fully developed, turbulent, two-dimensional channel flow where Rg 
was varied. His results indicate that, in general, the normalized turbulent 
quantities reduce as R§ is increased. Similar findings are obtained in the 
p rese nt measurements. As R§ is incr ease d by flow tripping, n' /\J^, v'/Uoof 
-n'v' -u'v'/u’v', and |u'v'j/q^ all show reduction, ccsnpared 
to the untripped flow where Rg is lower. 

The local skin-friction coefficients for both tripped and untripped bound- 
ary layers are found to be within 10 percent of the predictions made frcan 
Ludweig and Tillmann's empirical formula (eq. (6)), with the predicted values 
being lower. Although it was found that local distributions for mean or turbu- 
lent properties did not collapse for the two types of flow when the transverse 
coordinate was normalized to 6 , 6, or 0, similarity is observed for the wake 

velocity defect. Shown in figure 26 is the distribution of the wake-velocity 
defect with downstream distance normalized by the momentum thickness measured 
at the trailing edge. The similarity between the untripped and tripped flows 
is apparent. The asymptotic behavior of the wake-defect distribution also 
follows the expected (x/6)“®*^ variation. 


CONCLUDING REMARKS 

Measurements of mean and turbulent flow fields have been conducted near 
the trailing edge of a symmetrical NACA 63] -01 2 airfoil at zero angle of attack. 
The boundary-layer and near-wake characteristics were investigated for both 
naturally developed flew and artificially tripped flow. 

For the naturally developed flow, the general trends observed for the mean 
and turbulence profiles were found to be similar to those reported for a fully 
developed, turbulent boundary layer over a flat plate under zero pressure gra- 
dient. This suggests that the influences of the mild, adverse pressure gradient 
and the slight streamline curvature on the airfoil trailing-edge flew are weak. 
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Notable streamwise variations in turbulence properties were found in the airfoil 
flow field, however, which are not expected for a fully developed equilibriuni 
flow. 


The main effect of artificially tripping the boundary layer near the air- 
foil leading edge is to make both mean and turbulent profiles nearly independent 
of streamwise location. The changes observed in the mean- and turbulence- flow 
properties as a result of tripping are expected on the basis of the Re3^olds 
number based on the boundary-layer thickness. Reductions in skin-friction coef- 
ficient, turbulent intensities, Reynolds shear-stress, turbulent kinetic energy, 
and turbulence shear correlation coefficient occur when Reynolds number based 
on boundary-layer thickness is increased as a result of tripping. 

The changes in the turbulence properties observed for the untripped flow 
just downstream of the airfoil trailing edge are believed to be manifestations 
of the adjustment made by the coherent eddies in the boundary layer when 
encountering the trailing edge. This adjustment governs the production of 
trailing-edge noise. The fact that most of the observed changes in the turbu- 
lence properties for the untripped flow were absent in the tripped flow is 
believed to be attributable to the increase in fine-scale turbulence and a loss 
of correlation for the coherent eddies due to tripping. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
March 24, 1981 
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TABI.E 1 MEASUREMENT STATIONS 


Station no. 

x/c 

Flow region 

1 

0.9000 

Boundary layer 

3 

.9375 

Boundary layer 

5 

.9583 

Boundary layer 

6 

.9688 

Boundary layer 

7 

.9792 

Boundary layer 

8 

.9844 

Boundary layer 

9 

.9896 

Boundary layer 

ro 

.9948 

Boundary layer 

11 

.9979 

Boundary layer 

12 

1.0029 

Near wake 

13 

1 . 0052 

Near wake 

1 4 

1 .0063 

Near wake 

1 5 

1 . 01 56 

Near wake 

1 6 

1.0208 

Near wake 

1 8 

1 . 041 7 

Near wake 

20 

1.0625 

Near wake 

22 

1 .100 

Near wake 
















(a) Boundary-layer flow. 
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(b) Wake flow. 


Figure 10.- Mean velocity distribution; untripped. 
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(a) Boundary-layer flow. 
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(b) Wake flow. 

Figure 13.- Distribution of transverse turbulent intensity; untripped 
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(a) Boundary- layer flow. 



|b) Wake flow. 

Figure 14.- Distribution of Reynolds shear stress; untripped. 
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(a) Boundary-layer flow. 


sta., 12 13 14 15 16 18 20 22 


Y,mm 




(b) Wake flow. 

Figure 15.- Distribution of turbulent-kinetic-energy flux; untripped. 
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-u'v'/u'v' 

(a) Boundary-layer flow. 
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(b) Wake flow. 

Figure 16.- Distribution of turbulent-shear correlation coefficient; untripped 
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Figure 19.- Distribution of mean-flow integral properties; tripped flow. 
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(a) Boundary-layer flow. 
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(b) Wake flow. 

Figure 21.- Distribution of transverse turbulent intensity; tripped 





(b) Wake flow. 

Figure 22.- Distribution of Reynolds shear stress; tripped. 
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